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On Resistive Networks of Constant Power Devices 


John W. Simpson-Porco, Florian Dorfler and Francesco Bullo 


Abstract —This brief examines the behavior of DC circuits 
comprised of resistively interconnected constant power devices, 
as may arise in DC microgrids containing micro-sources and 
constant power loads. We derive a sufficient condition for all 
operating points of the circuit to lie in a desirable set, where 
the average nodal voltage level is high and nodal voltages are 
tightly clustered near one another. Our condition has the elegant 
physical interpretation that the ratio of resistive losses to total 
injected power should be small compared to a measure of 
network heterogeneity, as quantified by a ratio of conductance 
matrix eigenvalues. Perhaps surprisingly, the interplay between 
the circuit topology, branch conductances and the constant 
power devices implicitly defines a nominal voltage level for the 
circuit, despite the explicit absence of voltage-regulated nodes. 


I. Introduction 

The existence and uniqueness of operating points to non¬ 
linear resistive circuits is a classic topic in circuit theory, 
with nonlinearities usually entering in the form of voltage or 
current-controlled nonlinear conductances [1]. Many elegant 
approaches have been devised to study such nonlinear cir¬ 
cuits, from fixed point and global inverse function theorems 
to the topological concept of the degree of a mapping [2]. 
These approaches offer a binary answer to the question of op¬ 
erating point existence by inferring existence from continuity 
and monotonicity, or from sector boundedness conditions on 
the current-voltage characteristics. These results offer little 
guidance in estimating the locations of operating points in 
voltage-space, or in quantifying their behavior as a function 
of the circuit topology, branch conductances, and loads. 

In this brief we consider a modification of classic nonlinear 
resistive networks, with two distinguishing features. First, 
the nonlinearities considered herein arise not from nonlin¬ 
ear branch conductances internal to the circuit, but from 
externally-connected constant-power devices (CPDs). In con¬ 
trast to conductance models i = g{v) relating branch-wise 
voltage and current variables, an externally connected CPD 
constrains the relationship at the connection port between 
the port voltage V and the port current injection I. This 
hyperbolic constraint I = P/V results in a negative (resp. 
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positive) incremental conductance dl/dV for sources with 
P > 0 (resp. for loads with P < 0), and has been observed 
to lead to instability in power electronic [3], automotive 
[4], and power transmission systems [5]. Constant power 
loads are the most challenging component of the standard 
ZIP (constant-impedance/current/power) static load model, 
as their presence often renders analysis and design problems 
analytically intractable. Second, as a consequence of the 
CPDs regulating the power sinked or sourced through each 
port, the circuit lacks voltage-regulated nodes, and therefore 
lacks a nominal voltage level around which all nodal voltages 
cluster. For example, in power systems this nominal level 
is typically supplied by voltage-regulated buses such as 
generators and points of common coupling [5]-[7]. 

This work is motivated by the increasing deployment 
of microgrids. These small-footprint power systems offer 
reliability and flexibility by locally managing generation, 
storage, and load [8]. While microgrids can be directly 
connected to a utility, their major benefit comes from the 
ability to disconnect (or “island”) themselves and operate 
independently. Resistive circuits with CPDs arise in islanded 
DC and AC microgrids consisting of constant-power loads 
(CPLs) and micro-sources [9]-[ll]. As an example of a 
micro-source, photovoltaic panels are controlled for maxi¬ 
mum power point tracking [12], and appear to the network 
as constant sources of power. We refer the reader to [8]-[l 1] 
and the references therein for detailed modeling information 
on micro-sources and CPLs. Aside from this key application 
area and circuit-theoretic interest, the proposed analysis and 
its extensions may prove useful for novel forms of circuit 
reduction [13], synthesis [7], as a tool for reactive power 
flow analysis [6], [14], and in power transmission networks 
when generators reach their capability limits [5]. 

In Section |II] we develop the resistive circuit equations, 
and begin our analysis in Section by proposing a novels 
decomposition of the vectorized circuit equation. This de¬ 
composition can be thought of as separating the model into 
two equations: the first describes the resistive losses, while 
the second describes a complementary “lossless” flow of 
power. We present a simple and intuitive necessary condition 
for the existence of operating points. To build intuition for 
the general analysis which follows, in Section IV we study 
in detail the simplest CPD circuit consisting of two ports. We 
provide a necessary and sufficient condition for the existence 
of an operating point with positive voltages. In Section [V| we 
present our main result, partially generalizing the results of 
Section IV to arbitrary circuit topologies. We present suf¬ 
ficient parametric conditions for all circuit operating points 
to belong to an appropriately defined operating region, in 
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Fig. 1: A resistive n-port (n=3) with constant power devices. 
By convention, constant power loads satisfy 14/fc = Pk < 0- 


which the average nodal voltage is high and all voltages 
are tightly clustered near one another. A loose translation of 
the condition is the following: the ratio of resistive losses 
to total transmitted power should be small when compared 
to the (inverse) network heterogeneity, as quantified by the 
“eigenratio” [15] of conductance matrix. We regard our 
analysis as a first step in the theoretical understanding of the 
operating points of CPD circuits. While intuition suggests 
that the absence of voltage-regulated nodes will result in the 
circuit displaying a disorganized voltage profile, our main 
result shows that the voltage profiles of such circuits are 
quite uniform under appropriate conditions. 

The remainder of this section introduces some notation 
and preliminaries. Given a vector x G K", diag(a;) is the 
associated diagonal matrix. The n x n identity matrix is 
and (resp. 0„) is the n-dimensional vector of all ones 
(resp. all zeros). The set 1:^ = {x G K" : ija; = 0} is the 
subspace of all vectors in K" orthogonal to span(l„). For 
a vector x G K", ||a:||i = YJk=i Wl ||a ;||2 = {Y.'l=k 4)^^^ 
and ||a;||oo = maxfc \xk\- For a positive semidefinite sym¬ 
metric n X n matrix M, ||M ||2 = ^niax{M) (the largest 
eigenvalue), while ||M||oo = max^ \Mik\- 

II. Derivation of Network Equations 

In this paper we consider linear resistive circuits with 
nodes {0,..., n} for a total of n+l nodes. Nodes {1,..., n} 
are the ones of interest, which in our motivating example 
of a DC microgrid are locations where power will either be 
produced or consumed. We assume without loss of generality 
that these nodes form a connected network. For the short 
transmission lines we consider, the 0th node (ground) is 
electrically isolated from the others and taken as the datum. 
We may therefore consider the circuit as a “grounded n-port” 
where each port voltage is simply the node-to-datum voltage 
of the respective node [1]. After eliminating any passive 
internal nodes via Kron reduction [13], the input/output 
behavior of the n-port is described in the short-circuit 
admittance representation by [1] 

I = GV, (1) 


of the port current 1^ and port voltage I 4 to be a constant 
power Pk- The power injection Pk is positive for generation 
and negative for load. An n-port with CPDs is shown in 
Figure for a three-port circuit, where only one CPD has 
been shown. The port constraints due to CPDs read as 


P = diag(y)/, (2) 


where P = (Pi,..., P„). Substituting 0 into (|^, we arrive 
at our nonlinear network equations 


P = diag(E)GE. 


(3) 


For future use we collect some useful and well-known facts 
regarding the conductance matrix G [1]. 

Lemma II.l (Conductance Matrix) The n x n conduc¬ 
tance matrix G satisfies the following properties: 

(i) G = G^ is positive semidefinite; 

(ii) Gl„ = 0„; 

(hi) 0 = Ai(G)<A2(G)<---<A„(G). 


In contrast with nonsingular conductance matrix models 
derived from reduced node-edge incidence matrices [1] (i.e., 
incidence matrices in which the row associated to the datum 
has been removed), the conductance matrix G in our port 
model is singular due to the isolation of the datum node. 
This can be seen from Figure nodes {1,2,3} form a 
cutset, and hence the net current flow through this cutset 
— or equivalently, the sum over all port currents — must 
be zero. In vector notation, this reads as 1^1 = 0, or using 
the node equations ([T]), that l^GE = 0. This equality holds 
for all port voltages V if and only if Gl„ = 0 • 1„, which 
is Lemma II.l (ii). Equivalently, one may consider G as the 
conductance matrix of the sub-network containing all nodes 
other than the datum. 

The second eigenvalue A 2 (G) of the conductance matrix 
is called the algebraic connectivity [16], [17], and measures 
how strongly connected the circuit is. Similarly, the “eigen¬ 
ratio” A„(G)/A 2 (G) > 1 quantifies the heterogeneity of the 
circuit [15]. A large eigenratio A„(G)/A 2 (G) corresponds to 
a heterogeneous and/or weakly connected circuit, while an 
eigenratio close to unity corresponds to a highly symmetric 
and dense circuit with nearly uniform conductances. 

An operating point is any solution V of the circuit 
equations 0. In contrast to standard circuit problems, the 
circuit described by 0 has no voltage-regulated sources - 
all port voltages are free variables. The absence of voltage- 
regulated sources makes determining the operating points 
of (|^ challenging and non-standard. Intuitively, one might 
come to the conclusion that the inflexible power demands 
of the CPDs would conflict with one another, making equi¬ 
librium conditions impossible. In the next three sections we 
develop some basic theoretical results for the circuit model 
(|^ and find that this intuition fails. 


where I = (/i,...,/„) and V = (Ei,...,E„) are the 
vectors of port currents and voltages, and G is the n x n 
conductance matrix. Each port k of our n-port is interfaced 
with a constant power device, which constrains the product 


III. Decomposition of Circuit Equations 
The starting point for our analysis of ([^ is inspired by 
the following observations. Under normal conditions, circuits 
with voltage-regulated nodes (such as power transmission 
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Fig. 2: A two-port resistive CPD circuit. 


and distribution networks) possess “almost uniform” operat¬ 
ing points [6], where all nodal voltages are clustered around 
a nominal value Vq. We may express this as 

V = VoiK+x), (4) 

where Vb > 0 is a nominal voltage level and the vector x 
is a small, dimensionless deviation variable^ That is, the 
voltage profile is the sum of a uniform profile Vpln, plus a 
perturbation term described by x. While the solution-space of 
models such as ([^ is generally multi-valued and diverse [18], 
we nonetheless are interested in such “almost uniform” so¬ 
lutions, as these are the operating points relevant in practice 
[5]. Lemma [nT| shows that the conductance matrix naturally 
satisfies a related decomposition, since Gl„ = Inspired 
by these properties, we similarly decompose the vector of 
powers P = {Pi ,..., P„) in ([^ as 

P = ^l„ + Pi, (5) 

n 

where p|| G K and P± G 1:^. Such a decomposition is 
uniquely defined, as one may verify by noting that p|| = 
l^P = ^-L = nP, where 11 = 

In — -tn'^n '^^e projection matrix onto the subspace 1;^. 
Physically, p|| is the total power dissipated in the network 
(the net difference between sourced and sinked power), while 
Px can be roughly interpreted as the power which flows 
“losslessly” between nodes. Our first result shows that this 
change of variables allow us to decompose the power balance 
(|^ into two equations in orthogonal subspaces. 

Lemma III.l (Decomposition) The voltage vector V = 
Vo{ln + x) is an operating point of © if and only if 

P\\=yoX^Gx, (6a) 

Pl = Vq {Gx -b diag(a;)Ga;) — —1„ . (6b) 

Proof Substituting Q and (|^ into ([^ and using Lemma 
|II.l I Property (ii), one obtains 

Pj. + — In = Vq [In + diag(x)) Gx . (7) 

n 

Equation ( |^ is obtained by left-multiplying Q by 1J, while 
( |6b| l is obtained by left-multiplying Q by the projector H. 
The converse direction is immediate. □ 

A necessary condition for the existence of an operating 
point for (|^-(|6b]l is obtained by noting that the right-hand 

'without loss of generality, one may assume that x G 'if. 


side of ( |^ is nonnegative. This expresses the fact that the 
circuit is passive and must dissipate energy, that is, > 0. 

Proposition III.l (Necessary Condition) If p\\ < 0, then 
possesses no operating points. 

Before proceeding in Section M to a general analysis of the 
equations (|^-(|6bli, in Section ]rV] we build intuition for the 
system behavior by studying the simple but illustrative case 
of a two port circuit (Figure |^. 

IV. Example: Two Port Circuit 
The decomposed circuit equations (|6ali-(|6bli can be solved 
in closed form for the specific case of a two port CPD circuit. 


Theorem IV.l (Operating Point for Two Node Circuit) 

Consider the circuit © defined for two ports connected by 
a conductance g > 0. Then the circuit has a high-voltage 
operating point V = Vo(l„ -b x) with a high average 
voltage Vq and a small percentage deviation x if and only if 


P 1 + P 2 


P\\ 


e]o,i[. 


(8) 
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If this is condition holds, then the unique operating point is 


V = Vo 


-1 


(9) 


where 
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Proof Writing x = ■ [1; 1] and G = g • [1, -1; -1,1], 

the decomposition (|6a|)-(|6bl) reduces to 

P\\ = 45K)^a'Lx) 
ll-Pj_||i = 4gV[fa:inax, 

where p|| = -b P2 and ||Px||i = I]y=i \Pj-,i\ = \Pi “ 
P 2 I. The first equation is simply ( |6a| ), while the second is 
obtained by subtracting the two linearly dependent equations 
in ( |6bl l. Hence we calculate that Xmax = P||/||Fx||i and 
Vo = ±||Px||i/(2^pp|j') and we obtain (|^. □ 

Note from the condition that if Pi > 0, then P 2 < 0 
and vice versa. That is, one node must generate power if 
the other consumes power. Thus, ||Px||i = \Pi — F 2 I is the 
absolute sum of power injections/demands at the ports of the 
network, and the parametric condition ([^ admits the elegant 
physical interpretation that the resistive losses in the network 
P\\ — Pi + P 2 should be small compared to the gross power 
transferred through the networks ports (Figure |^. 

The dependence of both tCmax and Vq on the resistive 
losses p|| is intuitive: as p|| decreases, the deviation vector 
X becomes smaller and the average voltage level Vq rises. 
That is, the voltage profile becomes increasingly uniform. 
Conversely, as p|| increases, x increases in size linearly, while 
Vq decreases with the square root of the losses. We invite the 
reader to compare this result with standard load flow results 
for networks with fixed-voltage buses [19, Chapter 1]. Unlike 
the classic results in [19], the existence condition ([^ does 
not depend on the line conductance g. 















V. Sufficient Conditions for General Topologies 
We now present sufficient conditions on the circuit param¬ 
eters of a general n-port which guarantee that all operating 
points V = Vb(ln + t) of © have high average voltage 
levels Vo and small percentage deviations x, generalizing the 
“small losses relative to total power” intuition developed in 


Section IV Unlike Theorem IV. 1 we do not provide exact 
solutions for the circuit operating points, but only bounds 
for where in voltage-space operating points are located. The 
results are applicable to general n-ports of arbitrary topology. 

Theorem V.l (Operating Regions I) Consider the circuit 
equations 0 with the voltage and power decompositions 
@-(0 leading to the decomposed circuit equations (|6a|l- 
(6b|l. Assume that the circuit parameters satisfy 


A A 


P\\ 


■ g ]0 -[ (10) 

|P±I!2-P|| A2(G) ^ ^ 


and accordingly define 


V ■ — 

r m 1 n — 


1 - A 


A 

A 


P\\ 


A2(G) 
< 1 . 


> 0 , 


(11a) 


(11b) 


Then any operating point V = Vo(ln + x) of the circuit 
which exists satisfies Vq > Vmin and ||a;||oo < Tmax- 


Theorem | V. 1 1 can be understood methodologically as fol¬ 
lows. First, using the given data of the problem, compute A 
from the definition in ( [TOl l. If 0 < A < 1/2, then one may 
compute Vtnin and Tmax from (|1 la| l-( [TTb| ), and Theorem V. 1 
guarantees that all operating points V = Vb(ln + x) of the 
circuit will have a high average voltage Vo > Vmin with 
minimal deviations ||a;||oo < a^min- A graphical interpretation 


of Theorem V. 1 for a three-port circuit will be presented later 
in Section VI If A < 0 or A > 1/2, we can make no claims 


regarding the locations of operating points. 

Remark V.2 (Interpretation and Qualitative Behavior) 

The parametric condition ( [TOl i implie^that 2p||/||Pi|| 2 < 
A 2 (G)/A„(G), which may be interpreted as follows; the 
ratio of resistive losses to total power transmitted must be 
smaller than the spectral gap of the conductance matrix. 

It is useful to compare the sufficient condition ( [T0| ) to the 
exact two-port condition ([^ in Theorem IV. 1 by specializing 
( [T0| l for a two-port network. In this case it is known that 
A 2 (G) = A„(G) = g, and the condition ( fTO) ! reduces to 
F||/I1 -Pl||2 < 1/3. Thus, compared to ([^, the results of 
Theorem IV. II are more conservative due to the 2-norm and 
the factor of 1/3. Just as was observed in the two-port case 
of Theorem |I V. 1 1 the variables Tmax and 1/nin are closely 
related. In the regime of small losses P||/||F’_L|j 2 «l, we 
observe from ( [TOl l that A ~ « 1. It follows 

from dllab-dllbll that 


A2(G) 


v; 


/A2(G) ||P^||2 


P\\ A„(G) 


Y P|l A„(G)’ ~ ||Pt||2 A2(G) ■ 

^This follows from the fact that y <y/V — y) for y S [0,1[. 


Thus, the lower bound on the mean voltage level Vq scales 
with the resistive losses l/.^/pjJ, while the voltage percentage 
deviations x are bounded linearly with p\\. In this regime of 
low resistive losses, the mean voltage Vo is sensitive to any 
increase in network losses and will quickly decline, while the 
percentage deviations x will, at most, increase linearly. These 
observations lead us to conclude that, much like in standard 
power systems, an accurate balancing of power supply and 
power demand is crucial to the efficient and stable operation 
of CPD-dominated DC microgrids. □ 

Proof Since for any a; S 1;/ it holds that A 2 (G)||a;||| < 
cc^Gt < A„(G)||x|| 2 , from ( |6a| l we have that 


( 12 ) 


From ( |6bl l it holds that 

^P± = Gx + diag(x)GT - ■ 

Vg Vg n 

Taking norms on both sides and bounding, we obtain 
< A„(G)(||a ;||2 -f ||t|| 2 ) + ^ 2 ^F|I 


(13) 


"0 *^0 

where we have used the fact that ||diag(a ;)||2 = ||a;:||oo < 
l|a;||2- Now using the upper bound in ( [T^ and rearranging, 
([T3J becomes 


V/,> 
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P\\ 


A2(G) 

P\\ 

A2(G) 

p\\ 

A2(G) 

p\\ 

A2(G) 


1/2 
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We now calculate using ( |llb| i-( [TTal i and ( |T^ that 
1 


< hr 




P\\ 


< 


1 


P\\ 


A2(G) 


— Xn 


A2(G) Vtnin 

and hence ||a;||oo < Tmax as claimed. □ 

For completeness we present an analogous oo-norm con¬ 
dition, which similarly restricts the heterogeneity of the 
network by comparing the smallest branch conductance 
Pmin to twice the largest nodal degree ||G||oo- Depending 
on the particular topology and heterogeneity of the circuit 
under consideration. Theorem |V.3| may offer more or less 
conservative results when compared to Theorem | V. 1 1 


Theorem V.3 (Operating Regions II) Consider the circuit 
equations © with the voltage and power decompositions 
(§-(15) leading to the decomposed circuit equations (|6a|l- 
(6 b) . Assume that the circuit parameters satisfy 

IIGIIc 


A^ 


P|| _ 

|.^-L||c 30 P\\ t/min 


e ]o,i/2[, 


(14) 



























































Fig. 3: Operating point (red) and permitted operating region 
(blue) for three bus network. 


where Pmin > 0 is the smallest branch conductance, and 
accordingly define 


V ■ — 

F m 1 n — 


1 - A 


A 


Pll 

^min 


1 - A 


Then all operating points V = Vo(l„ + x) of the circuit 
satisfy Vq > V^in and ||x||oo < icmax- 


The proof of Theorem V.3 is similar to the proof of 


Theorem V.l and will not be reported here. 


VI. Case Study; Three Node Circuit 

To illustrate and test the results presented in Section |V] 
we consider the network of Figure [T] with parameters Pi — 
P 3 = -3kW, P 2 = 6 . 6 kW, gi 2 = 523 = IS, and gis = 0.5S. 
The (in this case, unique) nonlinear solution to (|^ is given 
by V = (201.4V, 205.2V, 223.6V), and hence Vb = 210.1V 
and X = (-0.04,-0.02,0.06). For these parameters, one 
calculates using ( [TOl i of Theorem | V. 1 1 that A = 0.12. The 
sufficient condition 0 < A < 1/2 is satisfied, and one readily 
calculates from ( |1 la| i and ( |1 lb| i that Vb > Vrnin = 122V and 
||a^||oo A aJmax = 0.14. The solution space is depicted in 
Figure which shows that the solution V lies inside the 
set defined in Theorem V.l as expected. While the bounds 


developed in Theorem V. 1 become increasingly conservative 


as the number of ports increases, they can be considered as a 
first step into understanding the fundamental physics of such 
circuits composed of constant power devices. 


VII. Conclusions 

We have examined the behavior of linear resistive circuits 
with constant power devices interfaced at the circuit ports, 
and have provided sufficient conditions on the network 
parameters for circuit operating points to have high average 
voltage levels and minimal differences in nodal voltage. A 
curious observation is that despite the absence of voltage- 
regulated nodes, the network lower bounds its own mean 
voltage level, as quantified by ( |1 la| i. 

While we have partially characterized the circuit operating 
points through bounds, it remains an open question whether 
the circuit equations 0 are exactly solvable for general 


topologies, and how the results presented change with full 
ZIP load models. As an outlook to control applications, the 
relative scaling of the voltage profile heterogeneity Xmax and 
the minimum average voltage Vmin reported in Remark [V.2| 
hint at strategies for voltage profile control in DC microgrids, 
in which renewable sources are optimized and dispatched to 
simultaneously adjust the coupled values P± and p| |. 
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